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Asymmetry of the electronic states in hole- and electron-doped cuprates: 
Exact diagonalization study of the t-t'-t"-J model 
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We systematically examine the asymmetry of the electronic states in the hole- and electron-doped 
cuprates by using a t-J model with the second-neighbor hopping t! and third one t" (the t-t'-t"- 
J model). Numerically exact diagonalization method is employed for a 20-site square lattice. We 
impose twisted boundary conditions (BC) instead of standard periodic BC. For static and dynamical 
correlation functions, averaging procedure over the twisted BC is used to reduce the finite-size effect. 
We find that antiferromagnetic spin correlation remains strong in electron doping in contrast to the 
case of hole doping, being similar to the case of the periodic BC. This leads to a remarkable electron- 
hole asymmetry in the dynamical spin structure factor and two-magnon Raman scattering. By 
changing the twist, the single-particle spectral function is obtained for all momenta in the Brillouin 
zone. Examining the spectral function in detail, we find a gap opening at around the k = (n, 0) 
region for 10% doping of holes (the carrier concentration x = 0.1), leading to a Fermi arc that is 
consistent with experiments, fn electron doping, however, a gap opens at around k = (n/2, n/2) 
and persists up to x = 0.2, being correlated with the strength of the antiferromagnetic correlation. 
We find that the magnitude of the gaps is sensitive to t' and t" . A pseudogap is also seen in the 
optical conductivity for electron doping, and its magnitude is found to be the same as that in the 
spectral function. We compare calculated quantities with corresponding experimental data, and 
discuss similarities and differences between them as well as their implications. 

PACS numbers: 74.25.Jb, 71.10.Fd, 74.25.Ha, 74.72.Jt 



I. INTRODUCTION 



High-temperature superconductivity in cuprates 
emerges with carrier doping into insulating cuprates 
classified as the Mott insulator. The carrier introduced 
into the cuprates is either an electron or a hole. Although 
the symmetry of superconducting order parameter is 
common with d wave in both cases ^££££& the phase 
diagrams exhibit asymmetric behaviors between the 
electron and hole carriers. The most prominent dif- 
ference appears in the antiferromagnetic (AF) region 
near the Mott insulator, where the AF order disappears 
with a small amount of carrier concentration (x ~ 3%) 
in a hole-doped cuprate L^-xSr^CuC^ (LSCO), while 
in an electron-doped cuprate Nd2-zCe :E Cu04 (NCCO) 
the AF order persists up to x = 0.15j2 A difference in 
magnetic properties is also seen in inelastic neutron scat- 
tering experiments: LSCO shows incommensurate spin 
structures for a wide range of x^ while in NCCO there 
is no incommensurate structures but commensurate ones 
are observed^ It is also an important difference that a 
spin-gap behavior observed in the underdoped region of 
hole-doped cuprates by the nuclear magnetic resonance 
experiments is not reported in electron-doped cuprates^ 

Differences in the electronic properties between the 
hole- and electron-doped cuprates are also observed in 
other experiments. The optical conductivity obtained 
from reflectivity measurements exhibits a pseudogap fea- 
ture at around 0.2 eV in the AF phase of NCCO^ but 
there is no such a feature in LSCO with the same carrier 
concentration. The occurrence of the pseudogap in the 
optical conductivity is correlated with the strong temper- 



ature dependence of the Hall coefficients and a metallic 
behavior in the c-axis resistivity. From angle-resolved 
photoemission (ARPES) experiments, it is clearly ob- 
served that hole carriers doped into the parent Mott insu- 
lators first enter into the k = (±7r/2, ±tt/2) points in the 
Brillouin zone and produce a Fermi arc^iA 2 ^ but elec- 
tron carriers are accommodated at around k = (±7r, 0) 
and (0, ±7r) and then the Fermi surface is formed in 
the superconducting region^ The doping dependence 
of core-level photoemission also shows different behav- 
iors of the chemical potential shift between NCCO and 
LSCO, 1 — These experimental data indicate the difference 
of the electronic states between hole- and electron-doped 
cuprates. 

The electronic states and magnetic proper- 
ties in the electron-doped cuprates have been 
theoretically examined by many groups, focus- 
ing on the comparison with the hole-doped 
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them, the present author has studied a t-J model with 
the second-neighbor hopping t' and third one t" (a 
t-t'-t"-J model) ^LiLiS By applying the numerically 
exact diagonalization technique based on the Lanczos 
algorithm to small clusters, the dynamical spin structure 
factor, optical conductivity, single-particle spectral 
function, and thermodynamic properties have been cal- 
culated, and it has been pointed out that the difference 
of AF correlations caused by the presence of t' and t" 
is a prime source of the contrasting behaviors in the 
electronic states between the hole- and electron-doped 
cuprates. 

In small clusters used in our previous works ^i&iS 
the momenta defined were discrete in the momentum 



space, since periodic boundary conditions (BC) were 
used. Thus, momentum-dependent quantities such as the 
single-particle spectral function suffer from the discrete- 
ness. Two-particle correlation functions also suffer from 
the finite-size effects under the periodic BC, because the 
two-particle operators are described as the convolution of 
the single-particle operators that are defined discretely in 
the momentum space. Therefore, it is necessary to intro- 
duce a method that overcomes such discreteness and to 
clarify whether such a method changes the conclusions 
derived from small cluster calculations under the peri- 
odic BC. In this paper, we introduce twisted BC for a 
20-site square t-t'-t"-J cluster. The introduction of the 
twist can make the momenta defined continuous in the 
Brillouin zone, and thus we overcome the difficulty of the 
discreteness in the spectral functioniM For two-particle 
correlation functions, we introduce an averaging proce- 
dure over the twisted BC. This procedure is known to re- 
duce finite-size effects. 3 - Therefore, results obtained un- 
der the twisted BC are expected to provide new informa- 
tion that has not been obtained under the periodic BC. 
The quantities examined in this paper are, in addition 
to the single-particle spectral function, several response 
functions in terms of spin and charge, i.e., the dynamical 
spin correlation, two-magnon Raman scattering, and the 
optical conductivity, together with several static correla- 
tions. 

Being consistent with the previous works under peri- 
odic BCfi£*i£ii& we find that AF spin correlation remains 
strong in electron doping, in contrast to the case of hole 
doping, in the presence of the second- and third-neighbor 
hoppings t' and t". This leads to a remarkable electron- 
hole asymmetry in the dynamical spin structure factor 
and two-magnon Raman scattering. The single-particle 
spectral function also show dramatic differences between 
hole and electron dopings. Along the nodal direction, 
i.e., the k = (0, 0)-(tt, it) direction, a quasiparticle band 
in the hole-doped system is gapless at the Fermi level as 
expected, while in electron doping the band is gapped 
up to the concentration x — 0.2. The gap is found to 
be correlated with the strength of AF correlation, indi- 
cating that the gap is magnetically driven. The spectral 
function around the antinodal region, i.e., k = (n,0) re- 
gion, shows contrasting behaviors: A gap appears in hole 
doping but not in electron doping, leading to a Fermi- 
arc behavior only in the hole-doped system. The gap is 
found to be sensitive to t 1 and t", as is the case of the 
nodal gap in electron doping. Such a Fermi arc behav- 
ior has not been detected under the periodic BC. In the 
optical conductivity, a pseudogap clearly appears in the 
electron-doped system after the averaging procedure over 
the twisted BC. The origin of the pseudogap is attributed 
to the strong AF correlation in the spin background. The 
gap is found to have the same magnitude as that in the 
spectral function along the nodal direction, indicating 
the same origin. In terms of pairing of carriers, we ex- 
amine the d-wave pairing correlation function. The pair- 
ing is found to be enhanced in the underdoped region of 



electron-doped system and also in the overdoped region 
of hole-doped one, being consistent with previous studies 
under the periodic and open BCi 36 i 37 Since the quantities 
examined have fewer finite-size effects as compared with 
those under the periodic BC, we can make a more pre- 
cise comparison between these results and experimental 
data. From the comparison, we discuss similarities and 
differences between them as well as their implications. 

This paper is organized as follows. In Sec. |Bj we in- 
troduce the t-t'-t"-J model and show outlines of the pro- 
cedure to calculate the single-particle spectral function 
as well as the correlation functions for a 20-site square 
lattice under the twisted BC. In Sec. Hill calculated re- 
sults of the doping dependence of magnetic properties 
such as the spin correlation functions and two-magnon 
Raman scattering are presented. Being consistent with 
experiments, AF correlation remains strong in electron 
doping. The single-particle spectral functions are shown 
in Sec. IIVI Asymmetric electronic states between hole 
and electron dopings are discussed, focusing on gaps that 
appear in different momentum spaces. Their implica- 
tions are discussed compared with experimental data of 
ARPES. In Sec. [V] the charge dynamics and pairing 
properties in the t-t'-t"-J model are discussed. The dop- 
ing dependence of the optical conductivity clearly shows 
asymmetric electronic excitations that are closely related 
to the single-particle properties. The d-wave pairing also 
shows remarkable asymmetric behaviors, and the role of 
the k = (it, 0) states for the pairing is discussed. The 
summary is given in Sec. IVII 



II. MODEL AND METHOD 

The Hamiltonian of a t-J model with the second- 
neighbor hopping t' and third one t" (a t-t'-t"-J model) 
reads 
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with d — x and y, 8' = x + y and x — y, and 6" = 2x 
and 2y, x and y being the unit vectors in the x and 
y directions, respectively. The operator c\, a = Ci ](7 (l — 
n i,-a) annihilates a localized particle with spin a at site 
i with the constraint of no double occupancy, and Si is 
the spin operator at site i. In the model, the difference 
between hole and electron doping is taken into account 
by the sign difference of the hopping parameters together 
with the difference of the localized particle^ For hole 
doping, the particle is an electron with t > 0, t' < 0, and 



t" > 0, while the particle is a hole with t < 0, if > 0, 
and t" < for electron doping. Although the ratios 
if /t and if'/t are material dependent^ we take in this 
study t'/t = -0.25 and if'/t = 0.12 for both the holc- 
and electron-doped cases, which is obtained by fitting 
the ARPES Fermi surface for NCCO with x = 0.15.22 
The \t\ is usually taken to be 0.35 eV£& We set J/\t\ = 
0.4. Hereafter, H = e = 1, and the distance between the 
nearest-neighbor sites in the two-dimensional lattice is 
set to be unity. 

In order to examine the single-particle spectral func- 
tion in the model, we use the exact diagonalization 
method for an iV-site square lattice with the translational 
vectors R a = Ix + my and R& = — mx + ly, being that 



l,m ^ and N = V 



If periodic BC are used for 



the lattice, the momentum ko for single particle is de- 
fined as k = 2ir(lni — rnn 2 )/A^x + 27t(to7j 1 — ln 2 )/Ny, 
n\ and n<i being integers that give discrete N points in 
the first Brillouin zone. Introducing BC with twist, we 
can define momenta continuously in the Brillouin zone^i 
This procedure gives smooth band structures even if we 
use finite-size lattices. The twist induces the condition 
that Ci+R^er = e^ a h,a and c i+ R b!CT = e'^Q^, with ar- 
bitrary phases <p a and 0b. Note that (f> a = (fn, = (it) 
corresponds to the periodic (antiperiodic) BC. Introduc- 
ing an arbitrary momentum k = k x x + n y y defining 
4>a(b) = K ' R-a(6)i the momentum for a given k reads 



k = ko + k . 



(2) 



In order for k to cover the full Brillouin zone, k needs 
to scan a square with the four corners that (k x ,k v ) — 
±tt/N(1 -m,l + m) and ±ir/N(l + m,-l + m). We note 
that imposing the twist is equivalent to transforming the 
operator c' ia Cj >(T in (JIJ into e 1K ^ Rj_R ^c! .Cj^ (T , Ri being 
the position of site i. The twist changes the hopping 
terms but not the exchange term in Q. 

For a given k, the single-particle spectral function 
A(k, uj) at zero temperature reads 



A(k,u>) = A_(k,uj)+A + (k,uj) , 



(3) 



with 
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where A + (A_) is the electron-addition (electron- 
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removal) spectral function. For hole doping, a^ a 
and Ck,<T for A + and A_, respectively, c ka (ck jC r) being 
the Fourier component of c[ (ci i(T ) with momentum k 
defined in (J2J and spin a. On the other hand, Ok.o- = Ck.<x 
and c k g for A + and A_ , respectively, for electron dop- 
ing. The ^o an d ^m represent the ground state with the 
energy E" and the final state with E^, respectively, for 
a given k. The chemical potential /i K is also dependent 
on k, which is defined as one half of the energy difference 



between the first ionization and affinity states of the sys- 
tem. In this study, we calculate A(k, u>) for a lattice 
with N = 20 (I = 4, m = 2) using a standard Lanczos 
technique with a Lorentzian broadening of 0.2|i|. The 
total number of k taken in the calculation is N K = 320; 
thereby, the Brillouin zone has 7r/40 meshes. 

In contrast to the single-particle spectral function, the 
momentum transfer in two-particle correlation functions 
such as the dynamical spin correlation function is re- 
stricted to discrete momenta defined by the A-site lattice 
even if we introduce the twisted BC. In order to evalu- 
ate various two-particle correlation functions under the 
twisted BC, we average the correlation functions over the 
N K points of k. This procedure is known to reduce finite- 
size effects^ 

The dynamical spin correlation function reads 

S ^^ = atEE K*m K\ *o)| 2 ^ - Kn + E%) , 
n m 

where S* is the Fourier component of the z component 
of the spin operator with momentum transfer q. A stan- 
dard Lanczos technique with a Lorentzian broadening of 
0.02 1 1 1 is used for each k in (JSJ. 

The two-magnon Raman scattering spectrum with the 
Big symmetry is given by 



w^EEi^^^i 2 ^-^ 



with the Raman operator for Bi g mode 

R = 2^(S i+x • S ; - S i+y • S ; ) 
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A standard Lanczos technique with a Lorentzian broad- 
ening of 0.1 1 1 1 is used for each k in (J7J). 

The real part of the optical conductivity under the 
electric field applied along the x direction is given by 





er(cj) = er sing (w) + <r Icg (uj) , 


(8) 


where the 


regular part <7 reg (<x>) is 




<7 re g(a;) = 


^E^EK^I^I^)! 2 ^-^ 


+K) 


with the x component of the current operator 
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The singular part <J s i ng (uj) in (|SJ) is related to the charge 
stiffness D, which is sometimes called the Drude weight, 
through 



; (w) = 2irDS(iu) 



(11) 



The D satisfies a sum rule 
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where the stress-tensor operator t xx is given by 

Txx = — t / y ( C iA-x..a C h°- + ^i-x,o- Ci .°' , 
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A standard Lanczos technique with a Lorentzian broad- 
ening of 0.1|i| is used for each k in calculating 0. The 
same broadening is used for the singular part ltTT|) . 

We also calculate static correlation functions as a func- 
tion of distance r. The spin correlation with staggered 
phase factors, which is a measure of the strength of AF 
correlation, is defined as 

Cspin(r) = t^ E ^ E <*? l S i+p ' S 'l *?> > ( 14 ) 



where the summation of p is taken to be bonds satisfing 
\p\ = r, and N r is the number of the bonds. The factor 
P r is 1 when the two sites are in the same sublattice, and 
— 1 otherwise. The charge correlation for doped carriers 

is 
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where the number operator ri\ is given by ri\ = 
^2 a c\ a Ci. a . In terms of superconductivity in the t-t'-t" -J 
model, the d-wave pairing correlation can be calculated, 
which is defined as 



w^iEmE^K/' 
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where Ai is the d a; 2_ J/ 2-wave singlet operator 

A ; = —f= y^(-l)°'(Ci, (r Ci +x ,- . + Ci.g-Ci-x.-o- 

V 2 cr 

III. MAGNETIC PROPERTIES 

In this section we show calculated results of doping 
dependences of the spin correlation functions and two- 
magnon Raman scattering for the N = 20 square lattice 
of the t-t'-t"- J model. 




FIG. 1: Staggered spin correlation C sp m(rO as a function of 
the two-spin distance r for an N = 20 t-t'-t"-J lattice, ob- 
tained by averaging over twisted BC. (a) Hole doping (t = 1, 
t' = -0.25, t" = 0.12, and J = 0.4) and (b) electron doping 
(t = -1, t' = 0.25, t" = -0.12, and J = 0.4). Solid squares, 
circles, upper triangles, and lower triangles are x — (Heisen- 
berg model), x = 0.1, x = 0.2, and x = 0.3, respectively. For 
comparison the t-J results are plotted with open symbols in 
(a). 



Figure Q] shows the staggered spin correlation C sp i n (r) 
in i|14|) for different carrier concentration x, where x = 
N c /N, N c being the number of carriers in the TV-site 
lattice. For the Heisenberg model (x = 0.0), the correla- 
tion is almost constant at large distance, indicating the 
presence of AF order. In hole doping, the AF correla- 
tion is suppressed and rapidly decays at large distance 
with increasing x. On the other hand, the AF correla- 
tion in electron doping is similar to that of the Heisen- 
berg model, although the magnitude decreases with dop- 
ing. This indicates the presence of AF order even in 
the electron doped systems up to x — 0.2. At x = 0.3 
the AF correlation almost disappears at large distance. 
We note that, although detailed studies on doping de- 
pendence of the magnetic correlation length for electron- 
doped NCCO have been reported) 40 i 41 the present results 
cannot be compared with them because of small system 
size. For a t-t'-t" -U model, where U represents the on- 
site Coulomb interaction, it has been reported^ 1 " that the 
temperature-dependent correlation length agrees with 
the experiments. 

The difference of AF correlation between hole and elec- 
tron dopings in the t-t'(-t")-J model has already been 
pointed out under periodic BCi 16 i 17 i 18 The present re- 
sults obtained by averaging over twisted BC confirm that 
such a difference in AF correlation is independent of BC 
and thus intrinsic to the model. The origin of the dif- 
ference comes from the sign difference of t' and t"m&l& 
Let us consider the Hilbert-space bases with Neel-type 
spin configuration in the spin background. The t' and 
t" do not change spin configuration of these bases be- 
cause of the same sublattice hoppings. This means that 
the self-energies of the bases are dependent on the values 
of t' and t". We can findiSii^ that the energies become 
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FIG. 2: Dynamical spin correlation function 5(q, w) for an 
TV = 20 t-t'-t"-J lattice, obtained by averaging over twisted 
BC. (a) Hole doping (t = 1, t' = -0.25, t" = 0.12, and J = 
0.4) and (b) electron doping (t = -1, t' = 0.25, t" = -0.12, 
and J = 0.4). Solid and dashed lines represent the data for 
x = 0.1 and x = 0.2, respectively. The dotted line at q = 
(tt, 7r) in (a) represents the data for the t-J model at x — 0.1. 
The momenta defined in the lattice are shown in (a). In (b) 
the edge of the spin-wave excitations in the Heisenberg model 
(x = 0) obtained by the linear-spin-wave theory is indicated 
by the downward arrow for each momentum. 



lower when t' > corresponding to electron doping. This 
stabilizes the Neel-type spin configuration; thereby, AF 
correlation remains strong in electron doping. For hole 
doping, on the other hand, the self-energies of the bases 
are increased and thus become comparable with other 
bases with different spin configurations. Such a mixture 
of various spin configurations gives rise to a similar effect 
caused by nearest-neighbor hopping of carriers. In fact, 
we find from Fig. ^a) that the spin correlation is similar 
to that of the t-J model. We note that the similarity be- 
tween the t-t'-t"-J and t-J model has not been observed 
in a cluster under the periodic BC, where the spin corre- 
lation decreases with increasing the magnitude of t' and 
t". 

The dynamical spin correlation functions JSJ in the 
N = 20 lattice are shown in Fig. [21 For x = 0.1 of 
electron doping, the excitation at q = (jr, tt) exhibits 
the minimum energy among the momenta defined in the 
lattice, and has the largest weight. Since the staggered 
spin correlation C sp i n (r) shown in Fig. ^ indicates the 
presence of long-range AF order, the finite-excitation en- 
ergy at q = (tt, 7r) can be due to the finite-size effect 
that inevitably causes a discrete energy separation be- 
tween the ground state and excited states. Away from 
(tt, 7t), the spectral weights are distributed at the higher- 
energy region, whose scale is comparable with the spin- 
wave excitation of the Heisenberg model whose lower- 
bound edges are denoted by the downward arrows. This 
again confirms that the ground state is the AF ordered 
state even in the presence of mobile carriers. With fur- 
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FIG. 3: Two-magnon Raman spectrum Ir(lu) with Bi g ge- 
ometry for an N = 20 t-t'-t"-J lattice, obtained by averag- 
ing over twisted BC. (a) Hole doping (t — 1, t' — —0.25, 
t" = 0.12, and J = 0.4) and (b) electron doping (t — —1, 
t' = 0.25, t" = -0.12, and J = 0.4). Solid and dashed lines 
represent the data for x — 0.1 and x — 0.2, respectively. The 
dotted lines represent the spectrum in the Heisenberg model 
(x = 0) obtained by averaging over the spectra for N =16, 
18, 20, and 26 lattices. 



ther doping of electrons (x — 0.2), the (7r,7r) spectrum 
loses its weight and the high-energy weights at other mo- 
menta shift to the lower-energy side, as expected from 
the reduction of AF correlation. This doping dependence 
is qualitatively consistent with recent inelastic neutron 
scattering measurements for an electron-doped material, 
Pr 2 _ x Ce K Cu0 4 .^ 

At a; = 0.1 of hole doping, the lowest-energy exci- 
tations are not at q = (tt, tt) but at (ir, 0) as shown 
in Fig. El^a), although the spectral weight is the high- 
est at (ir, it). For x = 0.2, the (tt, tt) weight decreases 
and the weight at (3tt/5, Att/5) shifts to lower energy. 
These behaviors indicate a tendency toward incommen- 
surate spin correlations reported in hole-doped materials 
such as LSCO (Ref. 8) and YBa 2 Cu 3 07-5. 43 At (tt, tt) for 
x = 0.1, we find, from the comparison between the solid 
and dotted lines, that the introduction of t' and t" shifts 
weights to lower-energy side. Since the equal-time spin 
correlations are almost the same, as seen in Fig.^ we can 
say that the effect of t' and t" on magnetic properties is 
the shift of AF spin fluctuation toward lower frequencies. 

Two-magnon Raman scattering also shows such con- 
trasting magnetic behaviors between hole- and electron- 
doped systems. Figure |21 exhibits the Bi g two-magnon 
Raman spectrum Ir(u>) in (JJJ for the N = 20 lat- 
tice. The dotted lines in the figure represent the Bi g 
two-magnon Raman spectrum of the Heisenberg model 
(x = 0), which is obtained by averaging over the spec- 
tra for N =16, 18, 20, and 26 lattices in order to reduce 
finite-size effects. The main peak due to two magnons ap- 
pears at around uo ~ 3J4^ For x — 0.1 of electron dop- 
ing, the two-magnon peak position remains unchanged 
with uj ~ 3 J as expected from the presence of AF order, 
though its weight decreases. In contrast, the two-magnon 



peak shifts to the lower-energy side in hole doping from 
u> ~ 3J to 2J at x = 0.1. At x — 0.2, there is no 
pronounced magnon peaks in both the hole and electron 
dopings. Such a contrasting behavior about the posi- 
tion of the two-magnon peak is consistent with recent 
experimental data, where in hole doping the peak posi- 
tion shifts to the lower-energy side^ but not in electron 
doping.— We note that the line shapes of the calculated 
spectra cannot be directly compared with those of ex- 
perimental ones because of excluding electronic Raman 
contributions. 



IV. SINGLE-PARTICLE SPECTRAL FUNCTION 

In this section, we present doping dependence of the 
single-particle spectral function A(k, uj) for both hole- 
and electron-doped t-t'-t"-J models, and discuss the sim- 
ilarity and difference between the calculated results and 
ARPES data. 



A. Half filling 

Let us start with the spectral function at half filling 
(x = 0). Figure 0] shows the weight map along the high- 
symmetry lines in the first Brillouin zone obtained by 
introducing the twist for the N = 20 t-t'-t"-J model as 
mentioned in Sec. \n\ To obtain continuous weights, we 
perform a smoothing procedure for each symmetry line 
that has 7r/40 meshes4£ The on-site Coulomb interaction 
that determines Mott-gap magnitude is set to be U/\t\ = 
4|i|/J=10. 

The top of the lower Hubbard band is located at 
k = (7r/2,7r/2). Quasiparticle energies at around (tt, 0) 
are lower than that of (tt/2, tt/2). The spectral weights 
at around (tt, 0) are suppressed in contrast to the case of 
the t-J model, where the quasiparticle energies at both 
k = (tt, 0) and (tt/2, tt/2) are almost degenerate and their 
weights are similar^ On the other hand, the quasipar- 
ticle at (tt, 0) in the upper Hubbard band is located at 
the bottom of the band. Therefore, the charge excitation 
with minimum energy is from k = (tt/2, tt/2) in the lower 
Hubbard band to (tt, 0) in the upper Hubbard band. In 
other words, the Mott gap in the two-dimensional insu- 
lating cuprates is an indirect gap, as previously pointed 
out by the exact diagonalization study of a Hubbard 
model with t' and t"4£ From such an indirect nature, 
we can expect that doped holes predominantly enter into 
the (tt/2, tt/2) region in heavy underdoping, while the 
electrons enter into the (tt, 0) region. 

We also see in Fig. 0] that the spectral weights at 
around the bottom of the upper Hubbard band are the 
largest among those at other regions. Since spectral 
weight at half filling is roughly proportional to the prod- 
uct of the weight of the Neel-type configuration in the 
Heisenberg ground state by that in the single-carrier fi- 
nal state, the large weights at around (tt, 0) imply that 
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FIG. 4: (Color online) Weight map of the spectral function 
for an TV = 20 t-t'-if'-J model at half filling along the high- 
symmetry lines. |£j = 1, t' /t = —0.25, t /t — 0.12, and 
J/\t\ — 0.4. Twisted BC are imposed on the lattice in calcu- 
lating the final states. For each BC a Lorentzian broadening 
of 0.2 |£| is used. The scale of the weight is shown in the bar 
at the right side of the panels. The on-site Coulomb inter- 
action that determines the Mott-gap magnitude is set to be 
U/\t\ = 4|tj/J = 10. The chemical potential is located at the 
zero energy. 



the Neel-type configuration is dominant in the final states 
near (tt, 0), and thus AF correlation is strong in the low- 
energy sectors of the single electron-doped system. Such 
a momentum-dependent feature comes from the presence 
of H and t"M 



B. Hole doping 

Figure [3] shows the weight map of the spectral func- 
tion for a two- hole doped system (x — 0.1) along the 
high-symmetry lines in the first Brillouin zone. The 
red lines in the figure represent a noninteracting tight- 
binding band with the same hopping amplitudes as t = 1 , 
t' /t = —0.25, and t" /t = 0.12, which is available for a 
guide of the band renormalization of the t-t'-t"-J model. 
From the weight map, we find that large spectral weights 
appear at around k = (tt/2, tt/2) just below and above 
the Fermi level. We also find that along the (0,0)-(tt,tt) 
direction the dispersion exhibits a slight downturn to- 
ward (0,0) at (tt/2, tt/2). On the other hand, at around 
(tt, 0) the spectra are located below the Fermi level with 
small weight and flat dispersion. These behaviors are 
consistent with the picture that doped holes predomi- 
nantly occupy the (tt/2, tt/2) region in underdoped sys- 
tem, which is expected from the dispersion at half filling 
as discussed above. 

The most interesting feature in Fig. [3] is a gapped be- 
havior near the Fermi level along the (tt, 0)-(7t, tt) direc- 
tion. This seems to correspond to the pseudogap ob- 
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FIG. 5: (Color online) Weight map of the spectral function 
for an N — 20 t-t'-t"-J model at the hole concentration x = 
1 - 18/20 = 0.1. t = 1, f = -0.25, t" = 0.12, and J = 
0.4. Twisted BC are imposed on the lattice. For each BC a 
Lorentzian broadening of 0.2t is used. The scale of the weight 
is shown in the bar at the right side of the panels. The red 
curves represent a noninteracting tight-binding band with the 
same hopping amplitudes. 
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FIG. 6: (Color online) Contour plot of A(k,uj = 0) for an 
N = 20 t-t'-t"-J model at the hole concentration x — 1 — 
18/20 = 0.1. t = 1, t' = -0.25, t" = 0.12, and J = 0.4. (a) 
Electron-removal spectrum A_(k, to = 0) and (b) electron- 
addition spectrum A+ (k, w = 0). The scale of the weight is 
shown in the bar at the right side of the panels. The red curves 
represents a noninteracting tight-binding Fermi surface with 
the same hopping amplitudes. 



served in ARPES experiments for hole-doped high-T c 
cuprates4& In order to see such a gap feature in more 
detail, we show in Fig. [5] the intensity map at the Fermi 
level (uj = 0) for both the electron-removal (A-) and 
electron-addition (A + ) spectra. Note that in finite-size 
lattices A_(k, uo = 0) is not equal to A + (k,uj = 0), 
in contrast with the case of the thermodynamic limit, 
and rather A_ (k, u> = 0) in this case is comparable with 
ARPES intensity near the Fermi level. In Fig.|Sfa), spec- 
tral weights are large near the nodal (0,0)-(ir,ir) direc- 
tion, and the weights decrease away from the nodal region 



FIG. 7: (Color online) Same as Fig.^Jbut for a t-J model. 



along the original noninteracting Fermi surface. This is 
similar to the so-called Fermi arc observed in the normal 
state of underdoped high-T c materials jASiiM 9 . We do not 
find pocket-like features in the electron-removal side at 
around (7r/2,7r/2), consistent with the experiments. 

We also note that the ridge of the spectra in Fig. HJa) 
reaches the momenta of (0.67T, 0) and (0,0. 6tt). Such 
a feature has experimentally observed in Na-doped 
Ca2Cu02Cl2r^ and theoretically demonstrated by using 
the method of equations of motion for the t-t'-J models 
In the present calculations, the presence of the weight at 
around (0.67r,0) and (0,0. 6tt) is interpreted as the effect 
of t' and t": The quasiparticle energy at around (ir/2, 0) 
is insensitive to t' and t" in contrast with that at (w, 0), 
as is expected from the tight-binding form. Thereby, 
the quasiparticle is kept close to the Fermi level at the 
(tt/2, 0) region. 

In the electron-addition side [Fig. EJb)], spectral 
weights predominantly spread along the (0, 7r)-(7r, 0) di- 
rection. In order to detect these weights, we need angle- 
resolved inverse photoemission with high resolution. 

The gap along the (7r,0)-(7T,7r) direction is sensitive 
to the magnitude of t' and t" . To see this, we show 
in Fig. [7| the spectral function for the t-J model with 
x = 0.1. We find a less clear gap observed along the 
(7r,0)-(7r,7r) direction, although the dispersion along the 
(0,0)-(7r,7r) direction is similar to that of the t-t'-t"- 
J model. Therefore, we can say that the long-range 
hoppings t' and t" are responsible for the formation of 
the gap and thus for a Fermi-arc behavior observed in 
ARPES experiments of hole-doped cuprates. 12,13,49 Re- 
cent ARPES experiments^ 1 " have shown that the flat 
band at around k = (tt, 0) is deeper in energy for 
Bi 2 _ z Pb z Sr 2 Cai_ K (Pr,Er) :E Cu20 8+5 (BSCCO) than for 
LSCO. Since t' and t" for BSCCO are known to be larger 
than that for LSCO^ the experimental data are consis- 
tent with the present picture that t' and t" predominantly 
control the pseudogap magnitude. 
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FIG. 8: Dependence of the gap energy -E gap on t' for an 
JV = 20 t-t'-t"-J model at the hole concentration x = 0.1. 
t — 1 and J = 0.4. -E gap is the minimum-energy difference 
between the electron-removal and electron-addition states at 
k = (tt, tt/4). The ratio of t" /t is kept at -0.12/0.25 = -0.48 
(solid squares) and (open squares). 



In order to examine this more quantitatively, we plot 
in Fig. |S] the t' dependence of the gap energy E geLp for 
two cases of t"/t' = -0.12/0.25 = -0.48 and t"/t' = 0. 
The £gap is defined as the minimum-energy difference be- 
tween the electron-removal and electron-addition states 
at k = (tt, 7r/4). At —t'/t = 0, the gap is finite but may 
come partly from the finite-size effect. With increasing 
—t'/t, -Egap increases for both the cases of t" /if = and 
—0.48. Compairing them, we also find that t" signifi- 
cantly contributes for the gap magnitude. Since if and 
if not only reduce quasiparticle energy at around (tt, 0) 
but also enhance low-energy AF fluctuation in the spin 
background, the magnitude of the gap seems to be re- 
lated to the strength of AF fluctuation. In fact, we find 
that the gap is also dependent on the value of J: The 
gap at -t'/t = 0.25 changes from 0.85* for J/t = 0.4 to 
0.51* for J/t = 0.2. 

From Fig.JSJwe see that gap opens not only along (tt, 0)- 
(tt,tt) but also along the (0,0)-(tt,tt) direction. In the 
latter, the gap appears in the electron-addition side (a; > 
0), and its magnitude is smaller than that at (tt, tt/4)iSS 
The gap value is also found to be dependent on if if Fig. [5j 
is compared with Fig. {7\ In order to detect this gap, we 
need angle-resolved inverse photoemission experiments. 

At x = 0.2 (four holes in the N = 20 lattice), the 
gap becomes less clear than that at x — 0.1, as shown in 
Fig. [5] Along the (0, 0)-(7T, tt) direction, the gap feature 
almost disappears, and a downturn of the dispersion at 
(tt/2, tt/2), seen at x = 0.1, vanishes completely. Further- 
more, a spectral distribution near the Fermi level along 
the (tt,0)-(tt,tt) direction becomes continuous, though 
the weight at around (tt, tt/4) is still small. We can say 
that overall behavior gradually approaches the noninter- 
acting band with increasing hole concentration. 
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FIG. 10: (Color online) Weight map of the spectral function 
for an N = 20 t-t'-t"-J model at the electron concentration 
X = 22/20 - 1 = 0.1. t = -1, if = 0.25, t" = -0.12, and 
J = 0.4. Twisted BC are imposed on the lattice. For each 
BC a Lorentzian broadening of Q.2\t\ is used. The scale of 
the weight is shown in the bar at the right side of the panels. 
The red curves represent a noninteracting tight-binding band 
with the same hopping amplitudes. 



C. Electron doping 

Figure ITUI shows A(k,ui) for a two-electron doped sys- 
tem (x = 0.1) of the N = 20 t-t'-t"-J model. The spectra 
are very different from those for hole doping. At around 
k = (tt, 0), an electron pocket is seen as expected from 
the spectral function at half filling (see Fig. Q}. Note 
that, in order to get such a clear pocket, the presence of 
AF order is necessary^ In fact, strong AF correlation 
remains in the lattice, as explained in Sec. lIIII Along the 
(0,0)-(tt,tt) direction, we find a clear gap in the disper- 
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FIG. 11: (Color online) Contour plot of j4(k,oj = 0) for an 
TV = 20 t-t-t -J model at the electron concentration x = 
22/20 - 1 = 0.1. t = -1, t' = 0.25, t" = -0.12, and J = 0.4. 
(a) Electron-removal spectrum A_(k, oj = 0) and (b) electron- 
addition spectrum A+(k, u) — 0). The red curves represents 
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FIG. 12: (Color online) Same as Fig^but x = 24/20-1 = 
0.2. 



The Fermi surface map is shown in Fig. 1111 where elec- 
tron pockets centered at (n, 0) and (0, it) are clearly seen. 
The presence of the pockets and the absence of weights 
along the nodal direction are consistent with ARPES 
datai^ for underdoped NCCO where AF long-rang or- 
der persists. We also note that recent spectral function 
calculations for a Hubbard model with long-range hop- 
pings display gap behaviors consistent with our results 
as long as U is largei 2 ^ 2 ?'?^? 1 '? 2 ^? 

The spectral function for a four-electron doped system 
(x = 0.2) is exhibited in Fig. H21 The (it, 0) electron 
pocket seen at x = 0.1 almost disappears, although the 
spectral intensity at around (n, 0) is still strong enough to 
show a remnant of the pocket. In contrast, the (n/2, ir/2) 
gap clearly remains but with smaller gap magnitude. 
With further doping, the spectra show dispersions sim- 
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FIG. 13: (Color online) Same as Fig.ESlbut x = 26/20-1 
0.3. 



ilar to a noninteracting system. Figure ^] exhibits the 
case of x = 0.3, where the gap at around (ir/2, tt/2) al- 
most disappears and the dispersion qualitatively follows 
the noninteracting band. The velocity at the Fermi level 
is almost one half of the noninteracting-band velocity, 
which is independent of the Fermi momentum. 

Let us discuss the origin of the gap along the (0, 0)- 
(ir, n) direction seen for x — 0.1 and 0.2. The gap is, of 
course, the consequence of the presence of t' and t", be- 
cause there is no gap at the Fermi level along the nodal 
direction in the t-J model (see Fig. [7J). In Fig. 1141 we 
show the t' dependence of the gap energy E gap denned at 
the noninteracting Fermi momentum k F along the nodal 
direction. Since the t-J model has no gap, the gap value 
of 0.2|/j| at —t'/t = is due to the finite-size effect. As 
is the case of hole doping (see Fig.0, we find that E gap 
increases with increasing —t'/t for both x — 0.1 and 0.2. 
From the comparison between the cases of t"/t' = and 
—0.48 for x — 0.1, we also find that t" significantly con- 
tributes to the formation of the gap. Since t' and t" 
enhance AF correlation as discussed in Sec. lIIII the mag- 
nitude of the gap in electron doping seems to be related 
to the strength of AF correlation. Reflecting the sensi- 
tivity to the AF correlation, when we increase J from 0.4 
to 0.6 keeping t"/t' = —0.48 for x = 0.1, E gap increases 
from 1.01|t| to 1.33|t|. 

Even for x = 0.2, the gap increases with increasing 
—t'/t. This implies that AF correlation is still strong 
enough to produce the gap. However, comparing with 
the experimental fact that at this concentration the d- 
wave superconductivity emerges at low temperature, the 
gap obtained in the present calculation may indicate an 
overestimate of the AF correlation. Thus, we may need 
to clarify the mechanism of the gap closing which makes 
the system a d-wave superconductor. We will discuss this 
in the following section. 
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FIG. 14: Dependence of the gap energy -E gap on t' for 
an electron-doped N — 20 t-t'-t"-J model, t = — 1 and 
J = 0.4. -Egap is the the minimum-energy difference between 
the electron-removal and electron-addition states at the non- 
interacting Fermi momentum k F along the (0, 0)-(7T, tt) direc- 
tion. For the electron concentration x — 0.1, the ratio of t" /t 
is kept at —0.12/0.25 = —0.48 (solid squares) and (open 
squares). For x = 0.2, t" ft = —0.48 (solid circles). 



V. CHARGE DYNAMICS AND PAIRING 
PROPERTIES 

In this section, we first discuss the difference of the op- 
tical conductivity between the hole- and electron-doped 
t-t'-t"-J models. Next, we show the d-wave pairing cor- 
relations together with charge correlations in the models. 



A. Optical conductivity 

FigurelT51shows the dependence of the optical conduc- 
tivity cr(u>) on the carrier concentration x for the TV = 20 
t-t'-t"-J lattice. At x = 0.1, there is a broad-peak struc- 
ture at around uj ~ t in addition to the Drude contri- 
bution centered at u = for both the hole and electron 
dopings. Such a broad-peak structure is known to be 
incoherent charge excitations accompanied by magnetic 
excitations^ This is physically characterized as an exci- 
tation from the AF ground state to an excited state where 
wrong spin bonds are created by the motion of carriers. 
As a result of the presence of the broad peak separated 
from the Drude contribution, a gap- like feature, i.e., a 
pseudogap, emerges at around uj ~ 0.5i. 

For electron doping, it has been discussed^ that the 
pseudogap is very sensitive not only to J but also t' and 
t": Increasing the absolute values of t'/t and t"/t, the gap 
increases in energy. Such a pseudogap feature in u{uj) 
has been clearly observed in electron-doped NCCO^ 
In Fig. ll5f a). the gap feature is also seen in hole-doped 
case at x — 0.1. Although the gap feature has not been 
clearly reported in the normal state of hole-doped LSCO, 
a broad peak can be seen at uj ~ 0.5 eV for x ^ 0.06. 53 




FIG. 15: Optical conductivity a(uj) for an N - 20 t-t'-t"- 
J model, (a) Hole doping (t = 1, t' = -0.25, t" = 0.12, 
and J = 0.4) and (b) electron doping (t = — 1, t' = 0.25, 
t" = —0.12, and J = 0.4). Dashed, solid, and dotted lines 
represent the carrier concentration of X — 0.1, 0.2, and 0.3, 
respectively. Delta functions are broadened by a Lorentzian 
with a width of 0.1|£|. Insets: The a; dependence of the Drude 
weight D as well as the integrated total weight K. 



The calculated broad-peak structure in hole doping prob- 
ably corresponds to the broad peak observed experimen- 
tally. 

At x = 0.2, a remarkable difference appears between 
hole and electron dopings in Fig. ED A pseudogap re- 
mains in electron doping accompanied by a peak at 
uj = 0.7|t|, while it disappears in hole doping. Since such 
a gap feature is related to magnetic excitations as dis- 
cussed above, the difference should reflect the difference 
of magnetic properties. In fact, AF correlation at x = 0.2 
behaves differently as already shown in Fig. E where the 
AF order is expected for electron doping, while the AF 
correlation length is very short (smaller than 2 lattice 
units) in hole doping. This clearly demonstrates the fact 
that charge dynamics is strongly influenced by AF spin 
correlation. At x = 0.3, <t(uj) in electron doping shows 
similar behaviors to the hole-doped case. This is reason- 
able because the concentration of X = 0.3 is enough to 
kill AF correlation. 

It is also interesting to compare the peak position in 
cr(u>) with the gap in the single-particle spectral function 
A(k, uj) discussed in Sec. IIVI In electron doping, the 
values of the gap -E gap at around k = (tt/2,tt/2) are 1.0£ 
and 0.74£ for x = 0.1 and 0.2, respectively (see Fig. ll4[l . 
These numbers almost agree to the peak positions in <j(oj) 
as shown in Fig. ll5f b). Such an agreement indicates that 
the pseudogap in u(ui) and the gap in A(k, uj) in electron 
doping have the same origin. Needless to say, AF spin 
correlation is the underlying cause of the gaps. 

In the inset of Fig. ED the Drude weight D as well as 
the integrated total weight K defined in (|12l) is plotted 
as a function of x. Both the weights increase with x. 
Compairing the hole and electron dopings, we find that 
D as well as K is larger in electron doping than in hole 
doping. Such an enhancement in electron doping, par- 
ticular for x ^ 0.2, is a consequence of the interplay be- 
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FIG. 16: Charge correlation Charge (?") as a function of the 
two-carrier distance r for an N = 20 t-t'-t"-J lattice, aver- 
aging over twisted BC. Solid squares: hole doping (t — 1, 
t' — —0.25, t" — 0.12, and J = 0.4). Open squares: electron 
doping (t = -1, t' = 0.25, t" = -0.12, and J = 0.4). 
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FIG. 17: t' dependence of d x 2_ y 2-wa,ve pairing correlation 



Cpair(r) at the largest distance r = \/T0 for an TV = 20 t- 
t'-t"-J lattice, t = 1 and t" = 0.12 for t' = -0.25, while 
t = -1 and i" = -0.12 for t' = 0.25. J = 0.4 for both cases. 
t' — corresponds to the t-J model. Solid squares, circles, 
and triangles are for x = 0.1, x — 0.2, and 0.3, respectively. 



tween the charge motion and spin background, where the 
AF spin background makes possible smooth sublattice- 
charge flows via t' and t"i& 



B. Charge correlation and <i-wave pairing 

In order to consider the pairing of carriers, we first 
show in Fig.^|the charge correlation function C c harge (?") 
that gives information on attraction between doped carri- 
ers. Compairing hole and electron dopings, we find that 
the nearest- neighbor attraction (r = 1) is stronger in 
electron doping than in hole doping, and vice versa for 
the long-distance correlation (r = vlO), which is irre- 
spective of carrier concentration x. The relatively strong 
nearest-neighbor attraction in electron doping is easily 
understood if we consider the fact that the AF order 
existing in electron doping is favorable for the pair for- 
mation gaining the exchange energy. For hole doping, 
such a force to attract two carriers is weak and thus the 
carriers spread over the whole system, leading to the en- 
hancement of long-range carrier correlation. 

Since the nearest-neighbor charge attraction is strong 
in electron doping, short-range pairing of two electrons is 
also expected to be strong as compared with hole doping. 
However, this does not automatically mean that super- 
conducting pair-pair correlation is strong. In Fig. 1171 
we show the d-wave pairing correlation C pa i r (r) at the 
largest distance (r — VlO) in the N — 20 lattice. We 
find that the d-wave pairing correlation is strongly en- 
hanced for electron doping (£' = 0.25) at x = 0.1, which 
is consistent with a density- matrix renormalization-group 
calculation. 36 With increasing x from 0.1 to 0.2, the cor- 
relation decreases rapidly. The enhancement of pairing 
correlation for t' — 0.25 predominantly comes from the 
enhancement of pairing itself as indicated in Fig. 1101 In 



the momentum space, the strong pairing originates from 
the large single-particle spectral weights near the Fermi 
level at k = (ir, 0), as shown in Fig. ^] This is easily 
understood if we express the pairing operator A, in l|17(l 
as 2/7V X)k( cos ^ - cosky)J2 <J (~ 1 )' T ^,crC-k,-cr- The d- 
wave operator has the largest amplitude at k = (w, 0) 
and thus the large single-particle occupation at this k 
seen in electron doping at x = 0.1 gives rise to large 
pairing interaction. 

The enhancement of d-wave pairing correlation in elec- 
tron doping is accompanied by an enhancement of AF 
correlation. We speculate that the AF correlation ex- 
ceeds the pairing correlation near half filling, i.e., the 
AF order overcomes the superconducting order. With 
increasing electron concentration, AF correlation weak- 
ens and finally pairing correlation may become dominant, 
resulting in a transition from AF to superconducting or- 
der as observed experimentally. Such a picture may have 
some connections with recent data for Pri_ a: LaCe a ;Cu04 
indicating the coexistence of AF order and superconduc- 
tivity in the vicinity of the transition^* 

An important point to notice further is that the AF- 
superconducting transition may be accompanied by a 
topology change of the Fermi surface from small to large 
ones. At the same time, the gap at the Fermi level along 
the nodal direction is expected to be closed in order for 
d x 2_ y 2-wELve superconductivity to be induced. However, 
in the present calculations, the critical electron concen- 
tration x, where the gap closes (x ~ 0.3), is higher than 
experimental values of the AF-superconducting transi- 
tion (x ~ 0.1—0.15) or a quantum phase transition 55 
(x = 0.165). The discrepancy may indicate the presence 
of additional effects that have not been included in the 
present t-t'-t"-J model: for instance (1) the x dependence 
of the parameter values, which has been incorporated 
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into the studies of a t-t'-t"-U model) 28 i 29 i 30 i 31 i 32 and (2) 
the effect of inhomegeneity, the presence of which has 
been reported in local probes such as muon-spin relax- 
ation^ and nuclear magnetic resonanceSiS experiments. 
In any case, we may need to clarify the origin of the dis- 
crepancy. This still remains as a future issue. 

In hole doping, the d-wave pairing correlation shows 
small value up to x — 0.2 (see t' = —0.25 in Fig. I17fl . 
However, it is enhanced at x — 0.3. This is consistent 
with a recent theoretical study^i where the enhancement 
is attributed to the increase of the occupation number at 
k = (-7T, 0) whose position approaches the Fermi level. 



VI. SUMMARY 

In this paper, we have examined the doping depen- 
dence of magnetic and electronic properties in the hole- 
and electron-doped cuprates by using the exact diago- 
nalization technique for the t-t'-t"-J model. In order 
to reduce finite-size effects in small-size lattice, twisted 
BC are introduced instead of standard periodic BC. For 
the calculation of correlation functions, we have averaged 
the results for various twisted boundary conditions. The 
single-particle spectral function has been obtained for all 
momenta in the Brillouin zone by changing the twist. 

We find that the fact that AF spin correlation re- 
mains strong in electron doping in contrast to the case of 
hole doping, which has been obtained under the periodic 
BCfiS*i2iiS does not change even if the averaging proce- 
dure over the twist is employed. This confirms asymmet- 
ric magnetic properties in the t-t'-t"-J model. This nec- 
essarily leads to a remarkable electron-hole asymmetry in 
the dynamical spin structure factor and two-magnon Ra- 
man scattering. The doping dependence of these quan- 
tities in electron doping is qualitatively consistent with 
recent experimental data^*^ indicating the justification 
for the use of the t-t'-t"-J model. 

Using the twisted BC, we have also uncovered dra- 
matic differences in the single-particle spectral function 
between hole and electron dopings. In hole doping, the 
quasiparticle band for x = 0.1 is gapless at the Fermi 
level along the nodal (0, 0)-(tt, 7t) direction, but a gapped 
behavior emerges near the antinodal region. The Fermi 
surface map shows a Fermi arc behavior, consistent with 
ARPES data in the underdoped cuprates. It is important 
to notice that the presence of t' and t" is essential to the 
Fermi arc. The gap near the antinodal region disappears 
at x — 0.2. In contrast, the gap appears near the nodal 
region in electron doping up to x — 0.2. The gap is found 
to be correlated with the strength of AF correlation, in- 
dicating that the gap is magnetically driven. In addition 
to the gap, an electron pocket is clearly seen at around 



(tt,0) for a; = 0.1. 

In the optical conductivity, a pseudogap feature clearly 
appears in the electron-doped system up to x = 0.2 under 
the averaging procedure over the twist. The origin of the 
gap is attributed to the strong AF correlation in the spin 
background. In fact, we find that the pseudogap has the 
same magnitude as that in the spectral function along 
the nodal direction, confirming the same origin. 

Comparing the calculated spectral function and opti- 
cal conductivity with recent experimental data^2*A^ we 
find that the presence of the electron pocket and pseu- 
dogap agrees with the experimental data but only for 
x < 0.15. In particular, the critical electron concentra- 
tion x where the gap closes is x ~ 0.3 for the present 
calculation, but this is higher than experimental values 
of the AF-superconducting transition (x ~ 0.1—0.15) or 
a quantum phase transition 55 (x = 0.165). This discrep- 
ancy may indicate the presence of additional effects that 
have not been included in the present t-t'-t"-J model. 
We may need to clarify the origin of the discrepancy but 
leave this as a future issue. 

In terms of pairing of carriers, the d-wave pairing cor- 
relation function is examined, and the pairing is found to 
be enhanced in the underdoped region of electron-doped 
system and also in the overdoped region of hole-doped 
one, consistent with previous studies under the periodic 
and open BC. 3 £*2I In electron doping, AF correlation is 
also enhanced in the same concentration. We thus specu- 
late that AF correlation exceeds pairing correlation near 
half filling, but with increasing electron concentration AF 
correlation weakens and finally pairing correlation may 
become dominant. Although the electronic states in the 
normal state of high-T c cuprates, including asymmetry 
between hole and electron doping, are found to be de- 
scribed well by the t-t'-t"-J model, the relation of the 
t-t'-t"-J model to the d-wave superconductivity in addi- 
tion to the competition between AF order and supercon- 
ductivity in electron doping remains to be resolved in the 
future. 
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